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^sj ' Abstract. The dynamics of interacting structured populations can be modeled by ^^ = Ai{x)xi where Xi € R" 

X = {xi, . . . ,Xk), and Ai(x) are matrices with non-negative off-diagonal entries. These models are permanent if there 

t*^^, exists a positive global attractor and are robustly permanent if they remain permanent following perturbations of 

Ai{x). Necessary and sufficient conditions for robust permanence are derived using dominant Lyapunov exponents 

Ai(/^) of the Ai{x) with respect to invariant measures /i. The necessary condition requires max^ Ai(/i) > for all ergodic 

measures with support in the boundary of the non-negative cone. The sufficient condition requires that the boundary 

Cn , admits a Morse decomposition such that max^ A^ (fj.) > for all invariant measures fj, supported by a component 

Cn ' of the Morse decomposition. When the Morse components are Axiom A, uniquely ergodic, or support all but one 

population, the necessary and sufficient conditions are equivalent. Applications to spatial ecology, epidemiology, and 

.. gene networks are given. 

Oh 
O . Appeared in Journal of Differential Equations, 248, 1955-1971 (2010) 

qh' 1. Introduction 

A fundamental issue in population biology is what are the minimal conditions to ensure the long-term survivorship 
for interacting populations whether they be viral particles, bio-chemicals, plants, or animals. When these conditions 
^ ! are met the interacting populations are said to persist or coexist. Since the pioneering work of Lotka and Volterra 
on competitive and predator-prey interactions, Thompson, Nicholson, and Bailey on host-parasite interactions, and 
Kermack and McKendrick on disease outbreaks, nonlinear difference and differential equations have been used to 
understand conditions for population persistence [1^5]. One particularly important form of persistence for determin- 
istic models is permanence or uniform persistence which corresponds to the existence of a global attractor bounded 
r^ • away from extinction of one or more species. When such an attractor exists, interacting populations are able to 
f^ I recover from "vigorous shake ups" of the population state. Permanence has been characterized from a topological 
perspective [6|48| and with average Lyapunov functions |9l4l2| for differential equation models of the form 

(1) —J- =Xif,{x) i = l,...,k 

where Xi is the abundance of population i, x — {xi, . . . , Xk), and fi{x) is the per-ca pita gr owth of population i. 

An important extension of the concept of permanence is robust permanence |13l4l6l | which requires that (H)) 
remain permanent after sufficiently small perturbations of the per-capita growth rates fi. The importance of robust 
permanence stems from the fact that all models are approximations to reality. Consequently, if nearby models (e.g. 
more realistic models) are not permanent despite the focal model being permanent, then one can draw few (if any!) 
conclusions about the persistence of the biological system being approximated by the model. One can view robust 
permanence as one crude form of structural stability for population models. For this perspective, it is not unexpected 
that there are permanent systems that cannot be approximated by robustly permanent systems [l7| . In [IJ, [l6| , 
criteria for robust permanence were developed with respect to the average per-capita growth rates / fi{x) diJ,{x) with 
respect to invariant probability measures fi supported on the boundary of the positive orthant (i.e. where one or 
more populations are extinct). Roughly, these criteria for robust permanence require that the average per-capita 
growth rate (i.e. invasion rate) is positive for some missing species. 

While ([1]) can account for many types of population interactions, it assumes that all individuals within a popu- 
lation are exactly the same. However, theoretical biologists have long recognized that different individuals within 

1 



JOSEF HOFBAUER AND SEBASTIAN J. SCHREIBER 



a population may be in different states (e.g. different sizes or ages, livin g in different spatial locations) and these 



differences can have important consequences for population dynamics |18l - l24| . To account for how these differences 
influence persistence of interacting populations, we develop criteria for robust permanence for the dynamics of k 
interacting structured populations. Roughly, these dynamics correspond to replacing Xi in ([l} by vectors and fi{x) 
by matrices. In section [2l we describe these models in greater detail and review some basic concepts from dynamical 
systems theory. In section |3l we introduce the structured analog of average per-capita growth rates and define the 
concept of robustly unsaturated invariant sets. Necessary and sufficient conditions for an invariant set to be robustly 
unsaturated are proven. For Axiom A or uniquely ergodic invariant sets, the necessary and sufRcient criteria are 
shown to be equivalent. In sectional we use Morse decompositions and the criteria for robustly unsaturated invariant 
sets to develop necessary and sufficient conditions for robust permanence. In sections [5] and [6l we provide several 
applications of our results to spatially structured ecological models, structured epidemiological models, and models 
of gene networks. 

2. Models and assumptions 

Let Xi denote the state of the i-th population that lies in the non- negative cone C^ of R"' . Define x = {xi, . . . ,Xk) € 
C to be the non- negative cone of R" where n = X]i=i ^i- Let C^ denote positive cone {xi & Ci : JJ- xj > 0} for 
population i and C+ = Yii Cf, the positive cone for the interacting populations. If x i— >■ Ai{x) is a map into rii x Ui 
matrices that describes the growth of population i, then the dynamics of the interacting populations are given by 

(2) ^=Aix)x, i = l,...,k 

Let x.t denote the solution of ([2]) with the initial condition x, and more generally the semiflow generated by ([2]). 
Before stating our assumption on ([2|), recall a few definitions from dynamical systems. Given sets / C R and K C C, 
let K.I — {x.t : t G I,x ^ K}. A set iiT C C is invariant if K.t = K for all i > 0, and forward invariant if K.t C K 
for all i > 0. The omega limit set of a set K CI C equals u:{K) = nt>o ^■\^^ °*^)- The alpha limit set of iiT C C is 
a{K) = nt<o ^-("OOj A- Given a forward invariant set K , i? C -ft' is an attractor for the semiflow x.t restricted to K 
provided there exists an open neighborhood U Q K of B such that uj{U) = B. The stable set of a compact invariant 
set K is defined by 

W'{K) = {x&C: uj{x) ^ and uj{x) C A'}. 

The semiflow generated by ([2]) is dissipative if there exists an attractor B with ^{B) = C. 
Throughout this paper, we make the following assumptions: 

Al: X I— >■ Ai{x) is continuous 

A2: Ai{x) is irreducible and has non-negative off-diagonal entries. 

A3: x.t is defined for alH > 

A4: x.t is dissipative with global attractor r{A). 

Al is a basic regularity assumption. The non-negativity of off-diagonal entries in A2 implies that there are no 
negative feedbacks between individuals of different states in population i. This assumption is meet for many types 
of structured models as discussed in sections [5] and [6l The irreducibility assumption of A2 is generically meet and 
implies that all individuals within a population can pass through all states. In Remark [1] we discuss how this 
irreducibility assumption can be relaxed. A3 ensures that the population dynamics are defined for all future time. 
A4 requires that population densities/abundances eventually are uniformly bounded, a condition that should be met 
for an biologically realistic model. 

We say that the semifiow of ^ is permanent if this semiflow is dissipative and there is a positive attractor B C C+ 
such that W^{B) C C+. We will denote this positive attractor by A(^). Permanence ensures that populations can 
recover from rare large perturbations and allows for a diversity of dynamical behaviors. 

For K C C and i5 > 0, define the S -neighborhood of K as 

Ns{K) ^ {x e C : \x ~ y\ < 5 ioT some y £ K}. 
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We define a S -perturbation of 0) to be a system of the form 



that satisfies assumptions A1-A4, Y{A) C Ns{T{A)), and \A[x) - A{x)\ < S for all x e Ns{r{A)). © is robustly 
permanent if there exist 5 > and e > such that all (5-perturbations of ([2]) are permanent, and d{A{A), C \ C+) > e 
for all (5-perturbations i.e., there is a common/ uniform region of repulsion around the boundary. 

3. Invasion of compact sets 

We begin by studying the linear skew product flows on T{A) x R"' defined by {x,y).t = {x.t, Bi{t,x)y) where 
Y(t) — Bi{t,x) is the solution to ^{t) = Ai{x.t)Y{t) with 1^(0) equal to the identity matrix. Our assumption 
that Ai is irreducible with non-negative off diagonal entries implies that Bi{t,x)Ci C C^ for all x and t > (see, 
e.g., I25|). A result of Ruelle [26|, Prop. 3. 2] implies that there exist continuous maps Ui,Vi : T{A) — > C^ with 
Ui(a;)| = \vi{x)\ = 1 such that 

• the line bundle Ei{x) spanned by Ui{x) is invariant i.e. Ei{x.t) — Bi{x,t)Ei{x) for all i > 0. 

• the vector bundle Fi{x) perpendicular to Vi(x) is invariant i.e. Fi{x.t) — Bi{x,t)Fi{x) for all t>0. 

• there exist constants a > and /3 > such that 

(3) \\B,{t,x)\F,{x)\\ < aexp{-pt)\\B,{t,x)\E,{x)\\ 

for all x e K and t > 0. 

3.1. Invasion rates. Given any x S r(A), define the invasion rate of species i at population state x as 

Xi{x) = limsup — In \\Bi{t, x)\\ 

t-¥oo t 

When Xi = 0, \i{x) provides an upper bound to the rate of growth of population i when introduced at infinitesimally 
small densities. Important properties of this invasion rate are summarized in the following two propositions. For 
instance, the second proposition implies if {x.t)i stays bounded away from zero (i.e. population i persists), then 
\i{x) — 0. Alternatively, if \i{x) < 0, then population i is doomed to extinction. 

Proposition 1. Xi{x) satisfies the following properties: 

• Xi{x) = limsupf^Qo jln\Bi{t,x)v\ for any v > in R"% and 

• Xi{x) = lml8up^^^ J J^{Ai{x.s)ui{x.s),Ui{x.s)) ds. 

Proof. To prove the first property, we first show that 

Xi{x) = limsup - In \Bi{t,x)ui{x)\. 

t^oo t 

To this end, let v G R"' be any non-zero vector. Since R"* = Ei{x) (B Fi{x), we can write v = aui{x) -f w with a G R 
and w G Fi{x). Equation ([3]) implies 

\B,{t,x)v\ < a\B.,{t,x)u^{x)\ + \B,{t,x)w\ 

< \Bi{t,x)ui{x)\ (a -I- aexp(— /3t)|w|) . 

Hence, 

limsup - In |i3i(i,a;)w| < limsup - In |i?i(i,x)Mi(x)| < Xi{x) 

t— >CX0 t t— S-CXD t 

for all non-zero vectors v e R"'. Since ||i?i(i,x)|| = sup\^i^i \Bi{t,x)v\, this inequality implies that Xi{x) = 
limsupj_^gc j\n\Bi{t,x)ui{x)\. Now let v G C^. Then, we can write v — aui{x) + w with a > and w G Fi. 
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Equation ([3]) implies 

\Bi{t,x)v\ > a\Bi{t,x)ui{x)\ ~ \Bi{t,x)w\ 

> \Bi{t,x)ui{x)\{a — aexp{—l3t)\w\) 

Since a > 0, this inequality implies that 

Xi{x) > limsup — In \Bi{t, x)v\ > limsup — In \Bi{t, x)ui{x)\ — Xi{x) 

t— >oo t t— >oo t 

To prove the second assertion, let bi(t,x) = ln\Bi(t,x)ui(x)\. Invariance oi Ui(x) implies that bi{t,x) is additive: 

bi(t + s,x) = hi\Bi{t + s,x)ui{x)\ —hi\Bi{s,x.t)Bi{t,x)ui{x)\ 

= In \Bi{s, x.t)ui{x.t)\\Bi{t, x)ui{x)\ 

= In \Bi{s,x.t)ui{x.t)\ +hi\Bi{t,x)ui{x)\ 

= bi{s,x.t) +bi{t,x) 

Additivity of bi{t,x) and the fact that bi{0,x) — implies 



d b,{t + s,x) -~bi{t,x) 

—b,{t,x) = hm 

at s^o s 

bi{s,x.t) d 



lim 



bi{s,x.t) 



s^o s ds 

{j^Bi{s, x.t)ui{x.t),Bi{s, x.t)ui{x.t)) 



s=0 



\Bi{s,x.t)ui{x.t)\'^ 
= {Ai{x.t)Ui{x.t),Ui{x.t)) 

The Fundamental Theorem of Calculus implies 

bi{t,x) = / {Ai{x.s)ui{x.s),Ui{x.s)) ds 
Ja 

and the second assertion follows. D 

Proposition 2. For the solutions of ([2]) we have: 

• if Xi > 0, then Xi{x) < 0. 

• if Xi{x) < 0, then \iTnt~^ooix.t)i — 

• i/limsupj^o^ K^^-Oil > 0' then Xi{x) = 0. 

Proof. First, assume that x satisfies Xi > 0. Since the semiflow is dissipative, there exists 7 > such that \x.t\ < 7 
for alH > 0. Proposition [1] and the definition of the skew product fiow imply 

Xi{x) = limsup — In |_Bj(i, 2;)a;j| = limsup — In |(a;.i)i| 

i— ^00 t i— fcjo t 

In 7 
< hm sup ~ 

t— >-oo t 

To prove the second assertion, assume that Xi{x) < 0. If Xi = 0, then the invariance of the faces of C imply that 
{x.t)i ~ for all i > 0. Alternatively, if Xi > 0, then Proposition [1] and the definition of the skew product fiow imply 

limsup- ln|(xi)i| = Xi{x)<0. 

t-^oo t 

The final assertion follows from the first two assertions. D 
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3.2. Invariant measures. We review some definitions from ergodic tlieory. Given a Borei probabiiity measure /i 
on C, tire support of fi, denoted supp(/i), is the smallest closed set whose complement has measure zero. A Borel 
probability measure /i is called invariant for ^ provided that / h{x) d^(x) = J h(x.t) d^{x) for all i > and for all 
bounded continuous functions /i : C — > R. An invariant measure fj, is called ergodic provided that /i(-B) = or 1 for 
any invariant Borel set B. 

For an invariant measure fi define the invasion rate of species i with respect to fi as 

^i{f^) = / {At{x)ui{x),Ui{x))dfi{x) 
Jc 

Proposition 3. Let /i be an invariant measure. Ai(/i) satisfies the following properties: 

• Xi{x) = limf_^oo J jQ{Ai{x.s)ui{x.s),Ui(x.s))ds exists for fi- almost every x. Moreover, if ^ is ergodic, Xi{x) = 
Xi{fJ.) for fi-almost every x. 

• if ^ is ergodic, then there exists I C {1, . . . , fc} such that /i(Jlie/ ^i) ~ ^ '^'^'^ ^i{l^) = foi^ '^^^ i ^ I. 

Proof. Let /i be an invariant measure. Define hi{x) = {Ai{x)ui{x),Ui(x)). The Birkhoff ergodic theorem and 
Proposition [1] imply 



1 /■* 

Xi{x) — lim - / hi{x.s)ds 

t^OO t Jq 



exists for /i-almost every x. Moreover, Xi{x) = Xi{fj,) /i-almost surely if /i is ergodic. 

Assume /i is ergodic. By ergodicity and invariance of the faces of C, there exists / C {!,..., fc} such that 
/^(Ilie/ ^t) — 1- Let K C Yiiei ^t be a compact set such that fi{K) > 0. The Poincare recurrence theorem and 
the Birkhoff ergodic theorem imply that there is x G Yiiei ^t such that x.tn G K for some i„ t oo and Xi{x) — A,;(/i) 
for all i. The third assertion of Proposition [5] implies Ai(/i) = for all i G /. D 

3.3. Robustly unsaturated sets. Let K C dC be a compact isolated invariant set for the flow x.t restricted to 
dC. K is unsaturated if W^{K) C dC and K is isolated for x.t. If K is not unsaturated, then K is saturated. K 
is robustly unsaturated for (0j if there exists 5 > Q such that the continuation of K for any (5-perturbation of ([2|) is 
unsaturated. 

Theorem 1. Let K be a compact isolated invariant set for x.t restricted to dC. If one of the the following equivalent 
conditions hold 

• for all invariant measures fj, supported by K 

(4) max Ai(^) > 

l<i<A,' 

• there exist pi > such that 

(5) Y. P»^»(/^) > 

l<i<k 

for all ergodic probability measures supported by K . 

then K is robustly unsaturated for (0). Alternatively, if x ^^ A{x) is twice continuously differentiate and K is 
robustly unsaturated, then Ol) holds for all ergodic measures fi supported by K . 

Remark 1. For some applications (e.g., the disease model considered in section lKT]) . it useful to relax the irreducibil- 
ity assumption A2. For instance, if there exists an open neighborhood U of K such that for each i, Ai(x) has a fixed 
off-diagonal sign pattern for all x € U, then Ai{x) can be decomposed into a finite number, say mi, of irreducible 
components. For each of these irreducible components, one can define X^p) = J^{Al{x)ul{x),ul{x)) d^{x) where 
Al(x) is the submatrix of Ai{x) corresponding to thej-th irreducible component of Ai{x) and ul{x) is the continuous 
invariant subbundle for the irreducible linear cocycle determined by A^x). If we define Xi{ii) = maxi<j<TO; X\{pL), 
then all of the assertions of Theorem [7] still hold. 
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Proof. To prove the first assertion of the theorem, we need the following lemma. We call an invariant measure /x for 
([2]) saturated if Ai(/i) < for all i. 

Lemma 1. LetSn be a non-negative sequence that converges to zero as n -^ cxd. If fin are saturated invariant measures 
for Sn -perturbations of (0), then the weak* limit points o/{/i„}^2 '•^ "^ non-empty set consisting of saturated invariant 
measures for (0). 

Proof. Let dn be a non-negative sequence that converges to zero. Let /i„ be saturated invariant measures for Sn 
perturbations, x — A"'(x)x, of i = A[x)x. By weak* compactness of Borel probability measures supported on 
Ni{r{A)), there exist weak* limit points of /x„. Let fi be such a weak* limit point. Since /i„ are supported by 
Ni{r{A)) for all n sufficiently large, /i is supported by Ni{r{A)). To verify that fi is an invariant measure for ([2]), 
let /i : C — J> R be a bounded continuous function. Let t > and e > be given. Choose n sufficiently large so 
that I /(_, h{x) d^(x) — j^ h(x) (i/i„(2:)| < e, | J^-, h{x.t) dfj,{x) — J^-, h{x.t) (i/i„(a;)| < e, and \h[x.t) — h{x"- .t)\ < e for all 
x e 7Vi(r(A)). Then | J^ h{x.t) - h{x) dn{x)\ is 

< I / h{x.t) dii{x) — I h{x.t) diin{x)\ + \ I h{x.t) dfin{x) — / h{x) dfi{x)\ 

Jc Jc Jc Jc 

< €+ I \h{x.t)-h{x''.t)\dnn{x) + \ I h{x"' .t) dnn{x) ~ I h{x)dfj.{x)\ 

Jc Jc Jc 



< 2e+| / h{x'^ .t) ~ h{x) dfj.n{x)\ + \ / h{x) dfj.n{x) - / h{x) dfj.{x)\ 

Jc Jc Jc 

< 3e 

where the last line follows from the invariance of /Lt„ for x'^.t. Since e > and t > are arbitrary, we have 
/(-; h{x.t) djiix) = /f, h{x) dfi{x) for alH > and all bounded continuous functions /i : C — )► R. It follows that /i is 
an invariant measure for Q. To see that /x is saturated, define hf{x) — (Af {x)u'^ (x) , uf {x)) where u"(x) spans the 



invariant one-dimensional bundle given by [26|, Prop. 3. 2] for ^ — A"{x)x. Since A2{x) — >■ Ai{x) uniformly for x S 



Ni{A{A)) as n — >• cx), [26|, Prop. 3. 2] implies that u'^{x) — >■ Ui{x) converges uniformly for x £ Ni{A{A)). Given e > 0, 
choose n sufficiently large so that \hf {x) — hi{x)\ < e for all x £ Ni{A{A)) and | J^-, h^{x) dii{x)—Jfj /i"(x) d^n{x)\ < e. 
Then 

Aj(Ai) < / \h,{x) ~h'l{x)\dfi{x) + \ / h'Kx) dfj.{x) - / h'^ (x) dfin{x)\ 
Jc Jc Jc 



+ / h'l{x)dfin{x) 
Jc 
< 2e 

Since e > is arbitrary, Ai(/i) < for all i and /i is saturated. D 

Using this lemma, we prove that if K is saturated, then there exists a saturated invariant measure ^ supported by 
K . Assume K is saturated. Work of Hofbauer and So [7,, Thm.2.1] implies that either there exists y G W^{K)r\C^ or 
K is not isolated for the unrestricted flow x.t. If there exists y £ W''{K)r\C^ , then for all i > define i^t — j Jq Sy.s ds 
where Sy,s denotes a Dirac measure based at the point y.s. Dissipativeness of @ and weak* compactness of the 
Borel probability measures supported by A(^) imply there exist tk ^>- oo such that the sequence vt^ converges in 
the weak* topology to a Borel probability measure /i with support in K. A standard argument implies that n is x.t 
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invariant. Define /li : C — > R by hi{x) — {Ai{x)ui{x),Ui{x)). Proposition [2] and weak* convergence imply that 

1 /■* 
> Ai(j/) = hmsup- / hi{y.s)ds 

t— !-00 t Jq 

= lim sup / hi dvt 

> lim sup / hidvt^, = Xi{ii) 

/c— ^oo J C 

for all i. Alternatively, suppose that K is not isolated for the semiflow. Then there exists a sequence of positive 
w-limit sets that accumulate on K. Let /x„ be a sequence of ergodic probability measures supported by these oj-limit 
sets. Propositions [1] and |3] imply that /i„ are saturated for all n. Applying Lemma [1] with (5„ = for all n implies 
that there exists a saturated invariant measure fi supported by K. Hence, we have shown that if K is saturated, then 
there exists a saturated invariant probability measure fi supported by K. Equivalently, Q holding for all invariant 
probability measures with support in K implies that K is unsaturated. 

Next, we show that if K is not robustly saturated, then there exists a saturated invariant measure supported by 
K. Indeed, suppose K is not robustly saturated. Then there exists a non-negative sequence S„ converging to zero 
and a sequence of saturated measures /i„ for ^„-perturbations of ^ with support in the continuation of K. Let /i 
be a weak* limit point of {fJ.n}^=i- Lemma [1] implies that /i is saturated. Moreover, since the continuation of K 
converges to iiT as (5„ —^ 0, /i is a saturated invariant measure for ([2]) supported by K. Hence, we have shown that 
if K is not robustly saturated, then there exists a saturated invariant measure supported by K. Equivalently, if ^ 
holds for all invariant measures with support in K, then K is robustly unsaturated. 

To see the equivalence of the conditions given by (JH) and ([S]), let A — {p ^ R^ : J^iPi — 1} ^^"^ notice that 

niin max Ai (/i) = min max ^ piXi{fi) 

i 

where the minimum is taken over invariant probability measures fi with support in K. The Minimax theorem (see, 
e-g-, 127|) implies that 

(6) min max Ai (/i) = max min > piXi{fi) 

i 

where the minimum is taken over invariant probability measures /i with support in K. Since min^ ^^piAi(/i) is 
attained at an ergodic probability measure with support in K, the equivalence of the conditions given by (|4]) and ([5]) 
is established. 

To prove the final assertion of the Theorem, assume x ^^ A{x) is twice continuously differentiable and there exists 
a saturated ergodic measure /i supported by K. Proposition [3] implies that there exists / C {!,..., fc} such that 
A^dliG/ ^t) — 1 ^iid Ai(^) — for all i £ /. Since K C 9C, {1, . . . , fc} \ / is non-empty. We will show that ([2]) is not 
robustly permanent by proving that for all (5 > there exists a (5-perturbation of ([2]) for which the continuation of K 
is saturated. Let (5 > be given. Choose V C NsiT{A)) to be a compact neighborhood of r(yl) such that V.t C vaW 
for ah i > 0. Let 77 > be such that W =: N^{r{A)) cV. Let p : C -> [0, 1] be a C°° function such that p{x) = 1 
for aU X e T{A) and p{x) = for aU x £ C\W^.Define i = (ii, . . . , ifc) by 

Ai{x) Hi e I 



^'■^^^ ~ \ A,{x) - ^Idp{x) a 1(^1. 

where Id denotes the identity matrix of appropriate dimension. Let x.t denote the semiflow of i = A(x)x. Since 
xl = x.t whenever x.[0,t] eC\W,it follows that vTt C intl/ for aU t > 0. Hence, T{A) cV <Z Ns{T{A))). We also 
have ll^(x) — ^(a;)|| < |. Therefore, x — A{x)x is a (5-perturbation of ([2]). By construction, x.t — x.t for all x £ K 
and t > 0. Hence, /j, is ergodic for the semiflow of i = A{x)x and Xi{fi) < — | for this semiflow and i £ {1, . . . , fc} \ /. 
Let C and O be the Lyapunov exponents and Oseledec regular points supported by /i for i = A{x)x (see, e.g.. 
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[28j for definitions). At each point x G O, the splitting of R" determines three subspaces: the stable subspace 
E'^i^x), the center subspace E'^{x) and the unstable subspace E'^{x). Proposition [3] and our choice of A imply that 
E'^^x) n C+ ^ 0. The Pesin stable manifold theorem |28J, Corollaries 3.17 and 3.18] implies that tangent to E'^{x), 
E'^{x) and E"{x) are locally cc.i- invariant families of C^ discs W|, W^ and W^ corresponding to the stable, center 
and unstable manifolds. The family of stable manifolds W^ is contained in W^{A,snpp{iJ.)). Since E^{x) fl C+ ^ 0, 
W| n C+ ^ for some x ^ K. Consequently, K is saturated for x — A{x)x. D 

The difference between the sufficient and necessary condition is that for the sufficient condition, (|3]) has to hold for 
all invariant measures supported by K, while for the necessary condition, (j3]) has to hold only for ergodic measures 
supported by K. Since the invariant measures lie in the convex hull of the ergodic measures, the sufficient condition 
can be more restrictive than the necessary condition. 

For three classes of invariant sets, the necessary and sufficient conditions coincide. Recall, a compact invariant set 
K for ^ is uniquely ergodic if K only supports one invariant measure. Recall a compact invariant set K is Axiom 
A if the flow of ^ restricted to K is transitive and hyperbolic (see, e.g., [29| for deflnitions). 

Corollary 1. Assume x i— )■ A{x) is twice continuously differentiable. If K C dC is a compact invariant set and 
either it is uniquely ergodic or Axiom A, then K is robustly unsaturated if and only if ^ holds for all ergodic 
measures fj, with support in K . 

Proof. If K is uniquely ergodic, then the assertion follows immediately. For an Axiom A invariant set, Sigmund [30, 
Thm.l] has proven that invariant measures supported by periodic orbits of K are dense in the set of invariant 
probability measures supported by K. In particular, any invariant measure can be approximated by an ergodic 
measure and, consequently, the result follows. D 

Another special case where the necessary and sufficient conditions coincide is when K supports all populations 



except one. For discrete-time models, this case was considered by [22 1. 

Corollary 2. Assume x i— >■ A{x) is twice continuously differentiable. Let K C rii>2 ^i ^^ '^ compact invariant set. 
Then K is robustly unsaturated if and only if Ai(/i) > for all ergodic measures supported by K . 

Proof. Since K C ni>2 ^t ^ Proposition [3] implies that Ai(/i) = for all ergodic measures /i supported by K and 
2 < i < k. The ergodic decomposition theorem implies that Ai(/i) = for all invariant measures /i supported by K 
and 2 < i < k. Therefore, for any invariant measure /i supported by K, ^ holds if and only if Ai (/i) > 0. The 
ergodic decomposition theorem implies Ai(/i) > for all invariant measures supported by K if and only if Ai(/i) > 
for all ergodic measures supported by K. D 

4. Morse decompositions and robust permanence 

To state the sufficient condition for robust permanence, we use a characterization of permanence due to Garay [8| 
and Hofbauer and So [7[ that involves Morse decompositions of the boundary flow. Conley 31 1 defined a collection 
of sets {Ml, . . . , Aim} to be a Morse decomposition for a compact invariant set K if 

• Ml, . . . , Mm are pairwise disjoint, compact isolated invariant sets for the flow of ([2]) restricted to K. 

• For each x € K there are integers r = r(x) < s = s{x) such that a(x) C Mr and uj{x) C M^. 

• If r{x) — s{x), then x G M^. 

Garay, Hofbauer and So [7|, |8| proved the following characterization of permanence. 

Theorem 2 (Garay, Hofbauer-So) . // {Mi, . . . , M.^} is a Morse decomposition for dC, then (0) is permanent if and 
only if each of the components Mi are unsaturated. 

Theorem [1] and [2] imply the following result: 

Theorem 3. // {Mi, . . . , M„j} is a Morse decomposition for dC and (CI) holds for each of the components of the 
Morse decomposition, then 1^ is robustly permanent. Conversely, if x i-^ ^{x) *s twice continuously differentiable 



ROBUST PERMANENCE FOR INTERACTING STRUCTURED POPULATIONS 9 

and (III) is violated by an ergodic measure supported by one of the components of the Morse decomposition, then ([2]) 
is not robustly permanent. 

Theorem [3] in conjunction with CoroUaries [T] and [5] yield a characterization of robust permanence for a class of 
structured models. 

Corollary 3. Let {Mi, . . . , Mm} is a Morse decomposition for dC. Assume x i— > A(x) is twice continuously differ- 
entiate and for each Morse component Mi one of the following assertions hold 

• Mi is Axiom A, 

• Mi is uniquely ergodic, or 

• there exists j G {1, . . . , fc} such that Mi C {a; G C : x; S> for I ^ j}. 

Then ((2]) is robustly permanent if and only if Q holds for all ergodic measures supported by UiMi. 

Our results also provide a structured analogue for characterizing totally permanent systems [32|. For an ergodic 
probability measure /i, define supp(/i) to be the subset / C {1, . . . , fc} such that /^dlie/ ^+) = 1- 

Corollary 4. The following statements are equivalent: 

• for all ergodic probability measures fi with support in dC, 

Ai(/i) > for all i G {1, . . . , fc} \ supp{p) 

• ([2]) and all of its subsystems are robustly permanent. 

Proof. Suppose the first statement holds. Then unsaturated condition given by ([S]) with p = (1, . . . , 1) and Theorem[3] 
implies that 1^ and all its subsystems are robustly permanent. The other direction follows from the necessary 
condition for robust permanence in Theorem [3] and the second assertion of Proposition [3] D 

To illustrate the broad applicability of Theorem |3l we develop applications to spatially structured models in 
section [5] and to a disease model and a gene network model in section [HI 

5. Applications: Patch models 

A fundamental application of our results are to spatially structured models, with k species dispersing between m 
patches J33l - l35| : 



(7) M^^^//(^.)+^df^^_e^^^ ^^l,...,k j = l,...,m 

where xj denotes the density of species i in patch j, x^ — {x\, . . . ,x|,) is the vector of species densities in patch j, 
// is the per-capita growth rate of species i in patch j, d^ > is the dispersal rate for species i from patch I into 
patch j, and ej is the emigration rate of species i out of patch j. Hence ef > ^; d/. We assume that the matrices 
{df )j^i are irreducible for each species i. It is then easy to write ([7]) in the form ^, with Xi — {xj, . . . , x™)' where ' 
denotes transpose. 

For a single species ([7]) generates a monotone fiow. Under mild assuraptions, e.g ., each f^ is decreasing and 
negative for large densities, there is a globally stable equilibrium, see, e.g., [36[ and [3 71 sect. 5.4]. This equilibrium 
will be positive if the invasion rate at the origin A(0), which is given by the leading eigenvalue (stability modulus) of 
the matrix {d^^ + /-'(O) — e-'), is positive. 

5.1. Two species. For two competing species, e.g., each f^ is decreasing with respect to both species and is negative 
at large densities, the dynamics are still monotone, and hence almost all orbits converge to an equilibrium. Robust 
permanence requires two conditions. First, Ai(0) > for both species in which case there are two single species 
equilibria Ei and E2. Second, the invasion rates A2(Ei) and Ai(E2) are positive. Whether there is a unique positive 
globally stable equilibrium for the robustly permanent system depends in a delicate way on the system parameters, see 
[S^l- For slowly dispersing populations (i.e. < e^ <^ 1 for all i,j), the robust permanence condition is particularly 
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straight forward to verify. Let Ef be the largest solution to x^//(x^) = i.e. the equilibrium attained by species 
i in patch j when there is no dispersal and no competitors. For slowly dispersing populations, Ei is close to the 
product {El, . . . ,-E™), and the matrix A^(Ei) is close to the diagonal matrix diag(/2(£'^))"L]^. Hence the invasion 
rate A2 (Ei ) is close to maxj f^ {E{ ). So we obtain: The system with two slowly dispersing competing species is 
robustly permanent if and only if 

(8) max/^(£;^)>0 and max/j'(^^) > 

j j 

Intuitively, robust permanence requires for each species there is at least one patch where it can persist. 

For predator-prey systems the situation is even easier. Since (under mild assumptions on the prey dynamics) 
the Morse decomposition of dC is given by two equilibria {0, Ei}, the spatial predator-prey system is robustly 
permanent, if the predator can invade at the prey equilibrium: A2(Ei) > 0. However, the global dynamics are likely 
to be more complicated. For instance, if in each patch there is a globally stable limit cycle (e.g., |39l44lj ) and dispersal 
is sufficiently slow (i.e. < e;^ <i; 1 for all i,j), our results about robust permanence and the theory of normally 
hyperbolic manifolds [42| imply that there is a positive m dimensional torus, which attracts almost all orbits in C, 
whenever maxj /2(£'i) > 0. 

5.2. Rock-paper-scissors dynamics. The Lotka-Volterra model of rock- paper-scissor dynamics is a simple model 
that is used as prototype for understanding intransitive ecological outcomes 43|, |4J| . Here, a simple spatial version 
of this dynamic is given by 



dx-[ 
~dt' 

dx^2 
IT 



= x{{l ~ x{ — P''x2 — a-'xg) 4- y. d{ Xi — e{x{ 

k 



M ^ xiil-l3^x{-a^xi-xi) + Y,di''^''--'-' 



^3 ^3^3 



k 



where a^ £ (Oi 1) and /S-' > 1 for all j. A more general version of this dynamic is presented in [4J]. Under the 
assumption that there is no cost to dispersal (i.e. e^ = ^j d/ for all i,j), the maximal invariant set in dC consists 
of the origin and a heteroclinic cycle connecting positive single species equilibria Ei, E2, and E3. For slowly 
dispersing populations (0 < e^ ^ 1 for all i,j), Ei is close to (xi, 2:2, 2:3) — (1, . . . , 1, 0, . . . ,0, 0, . . . , 0) and the 
invasion rates A2(Ei) and A3(Ei) are close to maxj 1 — a^ and maxj 1 — /?■'. For the other equilibria E^, the invasion 
rates of the missing species are also given by maxj 1 — a-' and maxj 1 — j3^ . Consider the Morse decomposition of 
dC given by and the heteroclinic cycle. Since Ai(0) = 1 > for i = 1, 2, 3, is robustly unsaturated. An algebraic 
computation reveals that Q holds for all invariant measures supported by the heteroclinic cycle (i.e. all convex 
combinations of the Dirac measures supported by the equilibria) if and only if 



Equivalently, 



max(l — a^) > — max(l — (3^) 
j j 



2 > min a-' + niin /?-' 
j j 



In particular, even if the heteroclinic cycle is attracting for each patch when the system is uncoupled (i.e. 2 < a^ + /3^ 
for all j), it can be repelling for the weakly coupled system. 
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5.3. Three species Lotka Volterra in spatially homogenous environments. Consider Lotka-Volterra dy- 
namics in a spatially homogenous environment: 

(9) ^=a;^(r, -^a,,a;i) + ^dfx^-e^a;^', z = l,...,/c; j = l,...,m 

We assume that dissipativity can be shown by a linear Liapunov hmction: There are q > such that ^^ Ciai^XiXs < 
holds for all x G R^ \ {0}. Then the weighted sum of densities across all patches ^^ Cix\ is decreasing for large 
densitities and hence (jH]) generates a dissipative semiflow with global attractor T. We further assume that for each 
patch J and each species i the immigration rate equals the emigration rate, i.e., 

(10) eHE^''- 

This guarantees that the set H = {.t G C : a;^ = x\ \/i,i, 1} of spatially homogenous states is forward invariant under 

®- 

Now consider the Lotka-Volterra dynamics without spatial structure 

(11) -T^ = Xiij-i-'^aisXs), i = l,...,/c 

s 

It has been shown [12i, ch. 16.1, 16.2] that for fc = 3, ([TT|) is robustly permanent (robust meaning here within the class 
of Lotka-Volterra systems) if and only if all equilibria on the boundary of Ri^ are unsaturated and whenever there is 
a heteroclinic cycle connecting the one species equilibria then this cycle is repelling. Equivalently the following four 
conditions hold. 

(i) there exists an interior equilibrium x (i.e. Ax = r with a; 3> 0); 
(ii) dct(-A) > 0; 

(iii) the 2 species subsystems are not bistable competition systems; 

(iv) if there is a heteroclinic cycle between the one species equilibria, say Ei — > E2 — ;• E3 — ^ Ei then the following 
inequality holds 

(12) A2(Ei)A3(E2)Ai(E3) > |A3(Ei)Ai(E2)A2(E3)| 



Furthermore, the proof in 12|, ch. 16.1, 16.2] shows that the boundary flow of (fTT|) has a simple Morse decomposition: 
if (iv) applies then there are two Morse sets, the origin is a repeller, and the heteroclinic cycle is the dual attractor 
(within QR'^); if (iv) does not apply then the finitely many boundary equilibria form a Morse decomposition of the 
boimdary flow. 

Theorem 4. Under the above assumptions, ([9]) for k — 3 is robustly permanent if and only if (i), (ii), (iii) and (iv) 
hold. 

Proof. Suppose ([9]) is robustly permanent. Since the nonspatial system (fTT|) is the restriction of ([9]) to the invariant 
subspace H, PT|) must be robustly permanent and hence (i) - (iv) hold. 

Conversely, suppose that (i) - (iv) hold. We show that the above Morse decomposition of the boundary flow of 
pT|) (now in H) is also a Morse decomposition for ([9]) restricted to the maximal invariant subset of dC which is 
contained in the spatial two species subsystems of ©. For this we show that in each 2 species system the orbits 
of ^ converge to H (For the one species subsystems this is obvious.) If a two species subsystem of pTjl has an 
internal equilibrium, say E12 then by assumption (ii), this is not a saddle but a sink, and hence aiia22 > ai2a2i- 
This implies that the 1-2 submatrix oi A is VL-stable [12|, ch. 15.3], and [45| implies that E12 is the global attractor 
for the spatially structured system (J9l) restricted to the first two species. 

If a two species subsystem of (|11[) has no internal equilibrium, then a one species equilibrium say Ei is the global 
attractor. Then Lemma [2] (see below) shows the corresponding result for ([9]). A similar proof applies to the case 
where the origin is the global attractor of a two species subsystem. 
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Finally observe that under the above assumptions ([9|)- (fT0l) . the invasion rates \i at spatially homogenous boundary 
equilibria are the same for (|9]) and ([TT]) . Since all boundary equilibria are unsaturated for (fTT|) . they are for (|9]) as 
well. The same applies in case (iv) to the heteroclinic cycle. Note that the condition ((T^ on the eigenvalues is 
equivalent to the fact that all invariant measures supported on the heteroclinic cycle are unsaturated, see [la . Ex. 

4.5]. 

Lemma 2. Suppose, in a two species Lotka-Volterra system (|lll) (k ~ 2) all interior orbits converge to the one 
species equilibrium Ei. Then all orbits of the spatial version ^ (k = 2) converge to the spatially homogeneous one 
species equilibrium Ei G H. 

Proof. If the system is of predator-prey or mutualistic type then Easting's [45j result applies again. So let us assume 
that the local interaction is given by the competition system 

dxi dx2 

(13) — = rixi[l - xi ~ ax2) -rr ^ r2X2(l ~ pxi - X2) 

with < a < 1 < /3. In this case, Kazuo Kishimoto (letter to JH, Oct 1987) has given the following argument. 
There exists a family of forward invariant rectangles contracting to Ei. This shows that all interior solutions of © 
(fc = 2) converge to Ei. D 

D 

Theorem 3] shows that a spatial network of identical patches is permanent, if the within patch dynamics is 
permanent, and three species Lotka-Volterra dynamics. It is not clear how to extend this result to more than three 
species. The crucial step in the three species case is that permanence precludes bistable two species subsystems. 
The spatial version of a bistable two species system allows plenty of stable equilibria outside H, see [33|. However, 



permanent four species Lotka Volterra systems may have two species bistable subsystems [12|, ch.16.4]. So there may 



be invariant sets on the boundary outside H which need to be unsaturated to make Q permanent. 



6. More Applications 

6.1. Disease dynamics with density-dependent demography. Gao and Hethcote [46| introduced a model 
of disease dynamics with density-dependent demography. Here we consider a variation of their model in which a 
population of size P has a per-capita birth rate b(P) and per-capita death rate d. To ensure that the population 
persists in the absence of the disease, we assume that there exists K > such that b{K) = d. To describe the disease 
dynamics, let S be the number of individuals susceptible to the disease, / the number of infected individuals, and 
R the number of individuals that have recovered from the disease. We assume that P = S + I + R. If the disease 
transmission is asymptotic, /3 is the contact rate between susceptible and infected individuals, 7 is the rate at which 
individuals recover from the disease, and m the mortality rate due to the disease, then the population dynamics are 
given by 



b{P)P -dP- ml 



dP 
'dt 

^ = /3(P-/-i?)^^-(7 + d + m)/ 
dt e + P 

^ = ,1-dR 
dt 
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where e > is a constant. It is useful to introduce a change of coordinates in which y ~ js and z = ^■. 

dP 

_ = {b{P)-d-my)P 

f = (p{l-y-z)^-l-Hl-y)-b{P)\y 



ly - ib{P) - "m-y) z 



dz 
'di 
Setting xi = P and X2 = {y, z)' where ' denotes transpose yields a structured model where Ai{x) = b ~ d — my and 



^^(^) _ (P{^ - V - ^)lTP -l-rn(\^y)- b{P) 



7 my — b{P)^ 

Since {O,{K,0,0)} is a Morse decomposition of the boundary dynamics, Theorem [3] and Remark [T] imply that this 
model is robustly permanent if and only if 6(0) > d and Pjrrn > 7 + »77, + (i. In particular, for e > sufficiently small, 
one requires that the basic reproductive number , is greater than one. 

Remarkably, this same criterion determines robust permanence of models with significantly more complicated 
boundary dynamics. For example, suppose the focal population is a predator species. If the prey has abundance 
N and exhibits logistic dynamics, f{N) is the per-capita predator consumption rate of the prey, and b{N) is the 
per-capita reproductive rate of the predator, then the dynamics become 

dN 

rN{l - N/K) - J{N)P 



dt 
dP 



(6(iV) -d-my)P 



^ = (p{l-y~z)-^--1-m{l-y)-b{N)\y 

dz 

— = 72/ - {b{N) - my) z 

where r > and X > 0, and f{N), b{N) satisfy /(O) — b{0) — 0. Define xi — N, X2 — P, x^ = (y, z)' . Let us assume 
that the disease-free community (i.e. i/ = 2; = 0) is permanent (i.e. b(K) > d) and let A be the global attractor in 
the interior of the P~N plane. Under these assumptions, a Morse decomposition of the boundary dynamics is given 
by {xi,X2,X3) — (0,0,0), (a:;i,a;2,a:3) — (i^, 0,0), and A. Since the equilibria (0,0,0) and (A', 0,0) are unsaturated, 
it remains to characterize A3(^) for any invariant measure fi supported by A. Proposition [3] and Remark [T] imply 

A2(m) = [ ib{xi) - d) dfi{x) =0 
Jc 

Hence J^ b{xi) d^{x) — d and 

A3(/i) — P d^{x) ~ J — m — d 

J e + X2 

Theorem [3] and Remark [T] imply that this structured model is robustly permanent if , > 1 and e > is 

sufficiently small. 

6.2. Gene networks. The repressilator is an oscillatory gene network based on three (or more generally an odd 
number of) transciptional repressors. Such a system has been genetically engineered in E-coli by Elowitz and Leibler 



47| . Mathematical models go back to |48l. l49|. T hev all involve concentrations of proteins and mRNAs. A modified 
model with auto-activation was suggested in [50| 

(14) x^ = (3{y^ - x^) 

(15) yi = aF{x)i - yi 
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with 

(16) F{x)i=Xig{xi,Xi^i) 



1+Xi+ pXi^l + KpXiXi^l 

Here Xi and yi are normalized concentrations of proteins and mRNAs belonging to gene i, p is the strength of 
repression and k, a cooperativity parameter. This is a structured system of type ([2]) 



\mj \ag{x^,x^-i) -1) \y^ 

In [Sfl Thm 3], it was shown that for n = 3, a > 1 and p > 1, system (J14l) has a heteroclinic cycle connecting the 
3 single gene equilibria (similar to the rock-scissors-paper dynamics of section 5.2). Moreover 

• the system is permanent, i.e., the boundary of R^ is repelling, if A + /x > 0, 

• this heteroclinic cycle is asymptotically stable ii X + p < 

where A > and p < are the invasion rates of the two missing genes, i.e., the leading eigenvalues of the 2x2 
matrix in p7)) evaluated at the appropriate single gene equilibria (for explicit expressions see [501, eqs. (127, 128)].) 
Theorem [3] implies the permanence condition implies the stronger conclusion of robust permanence. A similar 
result holds for n odd, but the heteroclinic cycle is then actually a heteroclinic network. 
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